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I. FORMULATION OF THE PROBLEM 



Let's consider 2D Ising model on the simplest regular lattice with the partition 



jET function 

O ■ 

QJ 1 ( M-lN-1 



ttf ! {o"=±l} I m=l n=l 



m = l, ...,M, n = l, ...,iV. (1.1) 



A pair of numbers (m, n) enumerate the vertex of the lattice which is at the inter- 
section of m-th column and n-th line of the lattice. First the partition function ( 11.11 ) 
was calculated by L. Onsager [H], and then by several another authors and methods 
(see, for example, and others references can be found in [5]). Taking into 

account that o 2 mn = 1, rewrite the partition function as 



1.2^ 



> 
in 

(N| M-lN-1 

z= z) n n [(ch#+sh# ) (ch0 + sh<9 • (T m , n °m,n+l) 

VO 1 { cr =± 1 l 111=1 n=1 

O . 

The right-hand side of the last equation is a polynomial in the variables {cr m ,n}, 
each variable o m ^ n being in the degree not higher than 4. As a result of summation 
in ( 11.21 ) all summands with odd degrees of each of these variables reduce to zero. 
Thus one should take into account only the summands which are proportional to 
a mn = a mn = °mn ~ 1 f° r an 171 an d n - ^ follows from ( 11.21 ) that the sum is equal 
to the sum of all closed contours (loops) on the lattice, generally with intersections 
and self-intersections. Let's call the elementary parts of loops connecting the nearest 
vertexes (say, (m, n) and (m + 1, n) or (m, n) and (m, n + 1)) by connections. 
Thus the length of the loop is equal to the number of its connections. Number of 
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connections coming to each vertex can be equal to 0, 2 or 4. Thus the partition 
function ( 11.21 ) can be expressed as 

Z = ( c he) 2 ^ M ^ N ^2 MN ^(th0)", (1.3) 

loops 

where v is the total number of connections of the loop. 
Here an another method of calculation of partition function (11.11 ) is proposed. 
Let's consider the Clifford algebra with MN generatives: 

lxl v + l y lx = ^ xyi x, y = 1, MN. (1.4) 

Matrices {7^} are Hermitian and I assume that their dimensions to be 2 MN / 2 x2 MN ^ 2 , 
that is the minimal possible one at given number of matrices fl. Algebra ( PI) 
implies that the trace of any odd product of 7-matrices is equal to zero, and 

tr 7 x 7y = 2 MN/2 5 xy , tijxlylzlv = 2 MN/2 (5 xy 5 zv - 5 xz 5 yv + 5 xv 5 yz ) , (1.5) 

and so on. By virtue of algebra ( 11.41 ) any even product of 7-matrices is reduced to 
the ±1 (versus the number of their permutations) or to the product of two by two 
different 7-matrices. According to ( 11.51 ) the trace of the product is equal to ±2 MN I 2 
in the first case and to zero in the last case. 

The theory of 7-matrices in the spaces of large dimensions as well as in Hilbert 
space can be found in |6| . 

Further we shall consider that the indices x, y and so on enumerate the nodes of 
the lattice and thus they are integer vectors in the plane: x = (m, n); it is assumed 
that the index m increases on the right while the index n increases up. Thus, the 
matrix 7 X = 7 m n is related to each node (m, n). We define also the elementary basic 
vectors ei = (1, 0) and e2 = (0, 1), so that x = mej + ne2. 

The following statement is valid: the partition function ( 11.11 ) or ( 11.21 ) can be 
rewritten also as 



Z = 2 MN / 2 (ch 26i M -^ N - l hi 

X [(A + ^7x7x+e 2 )(A + M7x7x+eJ] [(A + ^7 x+ei 7x+e 1+ e 2 )(A + ^7x+e 1 7x+2e 1 )] X • • • 
• • • X [(A + /i7x+e 2 7x+2e 2 ) (A + ^7x+e 2 7x+ei+e 2 

)] X 

(A + A i 7x+e 1 +e 2 7x+ei +2e 2 )(A + A i 7x+e 1 +e 2 7x+2e 1 +e 2 

)]x-..}. (1.6) 



x 



The expression in braces in ( 11.61 ) is some polynomial in 7-matrices. The statement 
is true if the traces of all monomials of 7-matrices in ( 11.61 ) are non-negative (the 
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positivity condition which is proved further) and 

ip ch6 . */; she . , 2th6 

A = cos — = , ii = sin — = . — > smw = - . - ^ xo . (1.7) 

2 x/ch^ 2 v^h2^ l + (th#) 2 v ; 



Note that the last expression (11.61) is obtained from (11.21) by substitutions a 



7 mjTl and X^{cr=±i} — > 2 MN I 2 tr. But in contrast to (11.21) the sequence of multipliers 
arrangement as well as their form are crucial in the case (11.61 ). Draw attention to 
the fact that in the product under the sign of trace in (11.61) at first the brackets are 
multiplied in series along the lines, and then the results of multiplications along the 
lines are multiplied in series along the column. 

To establish the coincidence of the right-hand sides of Eqs. (11.61 ) and (11.31) . it is 
sufficient to prove the positivity condition and to take into account the fact that 
under the sign of trace in (11.61) only that polynomials in 7-matrixes " survive" which 
contain each of 7-matrix in degrees 0, 2 or 4 (as well as in the case of variables a) 
and use designations (11.71) . 

Let's prove positivity condition with the help of mathematical induction method. 
During the process of calculations the numerical factors are ignored since only the 
question about the sign of the loop is significant here. 

We begin the calculation with verification of the fact that the matrix factor of 
elementary cell, corresponding to elementary loop binding elementary cell with ver- 
texes 

x, (x + ei), (x + ei + e 2 ), (x + e 2 ), (1.8) 
is equal to unit. Indeed, according to (11.61) this matrix factor is equal to 

(7x7x+e 2 ) (7x7x+ ei ) (7x+e 1 7x+e 1 +e 2 ) (7 x+e 2 7x+ei+e 2 

) = 1. (1.9) 

The further calculation is based on the possibility of presentation of loop matrix 
factor as a product of matrix factor of smaller loop, binding smaller number of 
elementary cells as compared with initial loop, and the matrix factors of elementary 
cells supplying smaller loop up to the initial one. Each of the equalities in Figs. 
1 a and 1 b represents graphically equality between loop matrix factor in the left- 
hand side and smaller loop matrix factor times elementary cells matrix factors. 
All loop matrix factors are calculated according to order in (11.61 ). In Fig. 1 solid 
lines describe loops while dotted lines divide loop interiors into elementary cells. 
Equations in Fig. 1 are established easily by direct testing. For example, for Figs. 
1 a 1 and 2 we have: 



1 I 

1 2 3 
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Fig. 1 a 



11 12 



14 



15 



16 



Fig. lb 

[(7174) (7172) (7273) (7376 ) (7475) (7576)] = 
= [(7174) (7172) (7275) (7475)] [(7275) (7273) (7376 ) (757e)] , 
[(7779) (777s) (7s7io) (79711) (710712) (711712)] = 

= [(7779) (777s) (787io) (79710)] [(79711) (79710) (710712) (711712)] , (1.10) 

and so on. Developing loops and their matrix factors following to the schemes in 
Fig. 1, one can realize matrix factors of any loops as the product of elementary cells 
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matrix factors which are equal to unit according to (11.91 ). Thus the equality of any 
loops matrix factors to unit is proved inductively. Hence the positivity condition is 
established. 

It is evident that each multiplier (A + /i7 x 7 x+e ) in the braces in ( 11.61 ) is the matrix 
of orthogonal rotation by the angle ip in spinor representation in MiV-dimensional 
Euclidean space in the plane marked by the pair of 7-matrices in the multiplier. 
Therefore, all product in the braces in ( 11.61 ) is the matrix of an orthogonal rotation 
in spinor representation in M iV-dimensional Euclidean space. It is known js| that 
this matrix can be expressed in the form 



U = exp Qw x , y 7 x 7 y 



-UJ- 



and 



'x,y 



, , 1 

— <Jx,y + ^x,y + <^x,z^z,y + 



1.11) 



1.12) 



The trace of the braces in ( 11.61 ). i.e. the trace of the matrix ZY, is expressed simply 
through the eigenvalues of real orthogonal matrix O x y . Let the set of numbers 

(pi, Pi, P2, ~P~2> • ■ ■ > Pmn/2, pMN/2) (1-13) 
form the complete set of eigenvalues of the matrix 0^ y . Then (see Appendix) 



MN/2 



tiu= Y[ 



k=i 



2ch 



In p k 



MN/2 

n 

k=l 




MN/2 



= n (v^+ 

k=l 



Pk = e 



i<t>k 



1.14) 



In statistical limit, when M, N — > 00, the problem of matrix C X; y diagonalization 
simplifies radically since the matrices C X;X+Z and cj XjX+z become dependent only 
on z at relatively small distances from boundary of the lattice. This property is 
known as translational invariance. Therefore the diagonalization of these matrices is 
performed by means of Fourier transformation, i.e. by passing to quasi-momentum 
representation. The following complete orthonormal set of functions on the lattice 
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is used for that purpose: 



P 



q = 



%(m) = 


1 

ipm i \ _ 

Vm ' w ~ 




1 

) = 1 

} Vn 


TTylVl — 


0\ ir(M A\ 
Z J 71 I 1V1 — L ±\ 


..,0, 


9-7T 

Z 71 


M 


M ' ' 


M' 


7T{N- 


2) tt(N - 4) 


.,0, 


2tt 


N 


N ' ' 


iv' "■' 


k) = tt k ( 


x) = V p (m)V q (n), 


k = 





Y *k(x)* k -(x) = «5 k k> «— Y *k(x)* k (x') = ,5, 



1.15) 



II. CALCULATION OF EIGENVALUES 



At first let's solve a particular problem: define for x = (m, n) the spinor rotation 
operator transforming only 7-matrices with indexes (x + m'ei) and (x + m'ei + 62), 
where m! = 0, ±1, . . . : 



'2.r 



U {n) = . . . [(A + /i7x7x+e 2 )(A + ^7x7x+ ei )] x 

X [(A + ^7 x+ei 7x+e 1+ e 2 )(A + ^7x+e 1 7x+2e 1 )] 

Here we have the ordered product of spinor rotation matrixes along only one line 
n. The corresponding orthogonal matrix is defined just as in (11.121) : 

According to definitions ( 11.61 ). ( 11.121 ) and ( 12.11 ) 

U = .U^^U^U^ 1 ^ . . . . O = Q( n - 1 ) Q( n ) Q( n + l ) . _ . 



(2.2) 



[2.3) 



In the right hand sides of Eqs. ( 12.31 ) the ordered products of matrixes corresponding 
to lines take place. 

Let's find the obvious expression for the matrix O^l with the help of ( 12.21 ). This 
calculation is based on the simple relations 



(A + ^7 x 7 x+e ) t 7x (A + ^7 x 7 x+e ) = (cos V>)7x + (sin 
(A + / U7 x 7 x+ e) t 7x+e (A + /i7 x 7 x+e ) = (cos ip)j x+e ~ (sin^)7 x , 



[2A) 

;2.5) 
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(A + ^7x7x+e) f 7z (A + /i7x7x+e) = 7z, z/x,x+e, (2.6) 

following directly from Eqs. ( 11.41 ) and ( 11.71 ). Either vector ex or vector e2 can be 
taken instead of a vector e in Eqs. ( 12.41 )- ( |2TB| ). 

Consider the first case when y = (m, n). Then in the left hand side of Eq. ( 12.21 ) 
all spinor rotation matrixes forming placed to the right from (A + /i7 y _ ei 7 y )^ 
and all spinor rotation matrixes forming placed to the left from the matrix 
(A + /i7 y _ei7y) are cancelled mutually since they do not catch on matrix 7 y . But 
the "facings" 

(A + /iTy-eJy) 1 ■ ■ • (A + ^7y- ei 7y) 

transform 7 y in accordance with Eq. ( 12.51 ) in which one must make the identifica- 
tions e = ei, (x + e) = y. So we obtain the matrix element 

0j£_ ei = -(sin^), y = (m, n). (2.7) 



y + e 2 y + ei+e 2 y + 2ei + e 2 




y-ei 



y + ei 



y + 2ei 



Fig. 2 



Further only Eq. ( 12.41 ) is used in which the identifications x = (y + m'ei), m! 



0, 1, . . . are made in series, and at each m! at first e = e2 and then e = ei are put. 
This process is pictured symbolically in Fig. 2 where the bold point denotes the 
position of the initial 7-matrix placed in "facings" in the left hand side of Eq. ( 12.21 ). 
and the arrows show where its image comes under the action of the described linear 
transformations. Every time the result is reduced to multiplication by (sin?/;) if the 
matrix is transfered up or to the right, and to multiplication by (cos?/;) if the matrix 
is left on the former place. The matrix transfered up further remains unchanged. 
Thus the following result is obtained: 

0^ y+mfei = [sin^cosVT'(cos</0 3 , C'e I+ e 2 = [sin ip COS nl ^ m ' 



m! = 0, 1, . . . , y = (m, n). 



(2.8) 
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(n) 

The matrix elements as C y + e2 z are calculated similarly. The process of this cal- 
culation is pictured in Fig. 3 where the bold point again denotes the position of 
7-matrix at the beginning of the process, and the arrows show the subsequent move- 
ments of the images of initial matrix. At the first step the relation (12.51) is used in 
which the identifications e = e2, x = y are made. The relation (12.41) is used at the 
all subsequent steps. As a result of simple calculation we obtain: 



y + e 2 y + ei + e 2 y + 2ei + e 2 




y + ei 



y + 2ei 



Fig. 3 



O 



in) 



v v+p, = (cos ib), O } ,, ,,-n , = —(sin?/;) 5 [sin?/; cos 

y+ e 2,y+e2 v tji y+e2, y+(m'+l )ei+e2 V tj i r 



o 



(») 

y+e 2 ,y+m 



, = — [sin?/;cos?/>] TO+1 , m! = 0, 1, . . . , y = (m, n). 



2.9) 



Now it is necessary to make partial diagonalization of the orthogonal matrix using 
its translational invariance. For that purpose we make the partial Fourier transfor- 
mation along the lines of 7-matrixes with the help of Eqs. (11.151 ): 



,(p) = ^ *p(m) 7m , n = t1(-p), 7n(p), tI'O" 



;2.io) 



in 



Let's pass to Fourier transformation along the lines in Eq. (12.21) . Using translational 
invariance we rewrite this equation in the form 



U^ lnl (p)U^ n) = 



^m, 1 ri\m+m' 1 n"^P\" i ) 



ln"(p)- 



(2.11) 



The matrix in square bracket in the right hand side of Eq. (12.111) is designated as 
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{ p f n ^ n „ and it is calculated with the help of Eqs. (E2j)-(E9j): 



O 



(») 



V 



1 








• 




1 








• 







a p 


bp 


• 


n' =n 





-b p 


Cp 


• 


n 1 =71+1 











1 • 





\ 



;2.i2) 



/ 



Here 



a. 



cos ip 



sin?/') 



1 — e ip [sin iJj cos 



[sin ip cos 
1 — e l P [sin -0 cos 



;2.i3) 



and diagonal elements a p and c p are settled on n-th and (n + l)-th places, corre- 
spondingly. The direct check shows that the matrix ( 12.121 ) is unitary. 
Thus, according to ( 12.31 ) 



Op; n' ,n" — 



= { . . . 0^- l) 0^0^ +l) ...} . (2.14) 



The product ( 12.141 ) of matrix of the kind ( 12.121 ) is calculated easily: 

0, for n' < -1 

O 



p; n,n+n' 



-bp), for n' = -1 . (2.15) 
Cpbpdp, for n' > 

To find the eigenvalues of orthogonal matrix ( 11.121 ) one must calculate the Fourier 
components of ( 12.151 ) along the column: 



Pp,q Op 7 q ^ ] Op; n,n+n'^ q(jl ) 



n c _l ^2 _ 



1 - e^b r 



'V -lp,q 

rjp^q = 1 - (e~ ip + e~ l<z ) (sin?/; cos?/') 
Equation ( 12.161 ) is obtained with the help of ( 12.131 ). 



;2.i6) 
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III. THE PARTITION FUNCTION 



It follows from Eqs. (11.141 ). (11.151 ) and (12.161 ) that free energy is proportional to 
the following integral 

( ^ MN f +7r f +7r 

ln< j II (v^+ V%2) f = J d Pj dgln( v ^+ v / ^) = 

MN r\ f + \ l y P , q + rj p , q \ MN p . 

dp I dqln ( — ?=== I = -ryr I dp dqm {V P ,q + rip,q) ■ 



47T 2 



47T 2 / 



The last equality in (13.11 ) is true because 

/+7T /*+7T f-\-TT /*+7T 

dp/ dqlnr] P:q = dp dqlnfj pq = 
7T i/— 7T «/ — 7T 1/ — 7T 

Prove the equalities (13.21 ). According to (12.161 ) 

= (1 - ac"*) (1 - ftc _<? ) , 
sin 2i(j 



(3.1; 



(3.2) 



« = ^(sm2^), A, = - • 
2 2 — (sin 2*^Je ip 



(3.3) 



and /? p e 4,7 = 1 only in the case (sin 2ip) = 1, p = g = 0. Therefore 

/+7T /"+7T 
dp / dqlnr] Piq = 
-7T «/ — 7T 

dgjy dp In (1 - ae-*) \ + / dp <| / d g In (l - P P e~ iq ) \ = 

-EM-"/ 



n=l 



'+7T P+TT P+TT P+TT 

dq / dpe^+ / dp/3; / - 1 

-7T J —TV J —IT J— 7T 



dge- m n=0. (3.4) 



The second equality in (13.21 ) is proved similarly. 
Eventually, we can write out the free energy with the help of formulas (11.61 ). (11.71 ) 
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{HHP, (EH) and ((O): 



F= -TlnZ = 

(ch 2#) 2 - (sh 26) (cos j9 + cos q) 

= -MiVrj In 2 - ln(l - x 2 )+ 

g— / dj3 / dgln [(1 + x 2 ) 2 — 2x(l — x 2 )(cosp + cosg)] j, x = th(9. 

(3.5) 



The last expression in ( 13.51 ) coincides with the free energy of 2D Ising model given 
in |3j. 

The temperature of phase transition is obtained directly from ( 13.11 ) and ( 12.161 ). 
Indeed, free energy has a peculiarity at the phase transition point. It is seen from 
( 13.11 ) that this takes place when {r) PA + fj P q) —> for a part of quasi- moment a, which 
occurs only if (see ( 12.161 )) 

sin?/; c = cos^ c = — p, p — > 0, q — > 0. (3-6) 

v2 

Thus with the help of ( 13.61 ) and ( 11.71 ) we find that critical temperature is found from 
the equation 

th fl e = i^fe = v /2-l. (3 . 7) 
sin ipc 
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APPENDIX A 



Here the formula ( 11.141 ) is proved. 



Let lvx\ Vxj , k = 1, . . . , MN/2, be the complete orthonormal set of eigenvec- 
tors of the matrix x ,y, so that the eigenvalue pk (~p~k) corresponds to the eigenvector 
(vxj • Further also the designation 

{v£\ 1$} = {v$ , a = l,...,MN 
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is used. We shall consider the introduced vectors as vector-columns and the upper 
indices T and ' denote the transposition and Hermitian conjugation of vectors and 
matrices. By definition 



(Al) 

(A2) 
(A3) 
(A4) 

(A5) 
(A6) 
(A7) 

possess all properties of fermion creation and annihilation operators. Indeed, in 
consequence of (11.41 ) and ( jA2f > 

[c fc , 4,]+ = Sk^, [cfc, Ck>}+ = [4> 4']+ = °- ( A8 ) 

According to the definitions (IA3I ) and (IA7j ) we have 

Tx^xa = (4, Ci, . . . , 4w/2' C AfiV/2) ■ ( A9 ) 

With the help of Eqs. (IA5I ). ( 1A8j) and (IA9I ) the quantity (IA6I ) is rewritten as 

MN/2 MN/2 



The given definitions imply the following formulas: 

/7 xa = < or U = (y®, v®, v®, . . .) , (/7^) a6 = <U 

(L^OT)^ = diag (p b p l5 p 2 , P 2! • ■ • ) = Dab- 
It is shown in [6] that 

(U^uU) ab = diag (lnpi, - In pi, lnp 2 , - lnp2, ■ • •) = 
Due to (IA3[ ) and (IA5[ ) we have 

^7x^x,y7y = \ (Tx^xa) A a6 (^ y 7y) • 

2 mn/2 x 2^/ 2 -matrixes 



4 = 7x4? > c fc = 7x4? 



7x^x, y 7y = \Y1 [ ln ^ { C k Ck - CkC l 
k=l 



k=l 



(lnp k )clc k - -lnp k 



. (A10) 



Equality (11.141 ) follows immediately from ( jAlOj ) since the calculation of the trace 
in terms of 7-matrixes is equivalent to the calculation of trace in terms of the 
corresponding fermionic operators ( 1A7[ ). 
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